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ABSTRACT 

In the last ten years 20 millisecond pulsars have been discovered in the globular cluster 
47 Tucanae. Hitherto, only three of these pulsars had published timing solutions. 
Here we improve upon these three and present 12 new solutions. These measurements 
can be used to determine a variety of physical properties of the pulsars and of the 
cluster. The 15 pulsars have positions determined with typical errors of only a few 
milliarcseconds and they are all located within 1.'2 of the cluster centre. Their spatial 
density within that region is consistent with a distribution of the type n(r) oc r^^, 
with a sudden cutoff outside 4 core radii. Two pulsars have a projected separation of 
only 0'.'12, and could be part of a triple system containing two observable pulsars. We 
have measured the proper motions of five of the pulsars: the weighted mean of these, 
Ha = (6.6 ± 1.9) masyr"^ and = (— 3.4 ± 0.6) masyr"^, is in agreement with the 
proper motion of 47 Tucanae based on Hipparcos satellite data. The period derivatives 
measured for many of the pulsars are dominated by the dynamical effects of the cluster 
gravitational field, and are used to constrain the surface mass density of the cluster. 
The pulsar accelerations inferred from the observed period derivatives are consistent 
with those predicted by a King model using accepted cluster parameters. We derive 
limits on intrinsic pulsar parameters: all the pulsars have characteristic ages greater 
than 170 Myr and have magnetic fields smaller than 2.4 x 10^ Gauss; their average 
characteristic age is greater than ^ 1 Gyr. We have also measured the rate of advance 
of periastron for the binary pulsar J0024-7204H, oj = (0.059 ± 0.012)° yr-\ implying 
a total system mass l.AtoiM-Q with 95% confidence. 

Key words: binaries: general — globular clusters: individual (47 Tucanae) — pulsars: 
general 



1 INTRODUCTION 



Among the globular clusters in the Galactic system, 47 Tu- 
canae (henceforth 47 Tuc) holds the record for the number of 
known puls ars: 20 to date ([Manchester et al. 199~ 



Robinson 



et al. 1995; Camilo et al. 200C). Camilo et al. estimated the 



total number of potentially detectable pulsars in this cluster 
to be at least 200. The known pulsar population of 47 Tuc is 
very different from the population in the Galactic disk: all of 
the pulsars have periods less than 8 ms, and 13 are members 



of binary systems. Radio images of the cluster (Fruchter & 



Goss 2000; McConnell & Abies 2000) show three scintillat 



ing poi nt sources, the positions of which coincide with the 



brighte it known pulsars in the cluster: 47 Tuc C, D and J. 
It is therefore improbable that any bright pulsars remain to 
be discovered, even with the severe selection effects against 



the detection of binaries with very short orbital periods and 
pulsars with very short rotational periods noted by Camilo 
et al. (2000). 

In addition to this set of pulsars, there is a collection 
of other exotic objects near the core of 47 Tuc: at least 9 
X-ray sources (Hasinger, Johnston & Verbunt 1994; Ver- 
bunt fc Hasinger 1998) an d more than 20 blue stragglers 
([Guhathakurta et al. 1992). Evidence has also been found 



for large numbers of cataclysmic variables in the core (Ed- 
monds et al. 2001, in preparation). The high stellar density 
necessary for the formation of these obje cts may increase 
the rate of ex change of stars in binaries (Hut, Murphy & 



Verbunt 1991), leading to the formation of low- mass X-ray 



binaries and millisecond pulsars, a process addressed by Ra- 
sio, Pfahl & Rappaport (2000) for the specific case of 47 Tuc. 
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Table 1. Parameters of the globular cluster 47 Tuc. 



Parameter Value 

Centre R.A., a^r tuc (J2000) 00'^24™05';29 ± 0';28 

Centre Decl., 547 tuc {J2000) -72°04'52'.'3 ± I'.'S 

Distance, D 5.0±0.4kpc 

Age 10.0 ± 0.4 Gyr 

Total mass (1.07 ± 0.04) X 10^ Mq 

Tidal radius 40' (58 pc) 

Core radius, 9^ (rc) 23^1 ± 1^7 (0.6 pc) 

Escape velocity 58kms~^ 

Central density ~lx IO^Mqpc"^ 

Central line-of-siglit velocity dispersion, i'z(O) 11.6 ± 1.4kms~^ 

Proper motion in R.A., fia 7.0 ± l.Omasyr"^ 

Proper motion in Decl., fig —5.3 ± l.Omasyr"^ 



Reference 



( De Marchi et al. 199' 



( Dc Marchi et al. 199' 



( Giratton et al . 1997) 



( Vleylan 1989 j 



( Da Costa 1979|' 



( Howell. Guiiatliakurta fc Gilliland 2000| ) 



( Webbink 1985 ) 



fPrvor Mevlan 1993 



(Vlevlaii Mavor 1 9Sf ), 



f Ddenkirchen et al. 1997) 



( Jdenkirchen et al. 1997) 



Table ^ summarises the basic parameters of 47 Tuc 
which we use to interpret our results. This compilation has 
five important updates on the information tabulated by 
Camilo et al. (2000). First, the position of the centre of 
47 Tuc is taken to be point "D" of plate 1 of De Marchi et 
al. (1996). We derive the absolute position of "D" by com- 
paring its position with that of point "G" on the same plate, 
whose absolute position was obtained by Guhathak urta et 
al. (1992). Second, recent Hipparcos data (Reid 1998) yields 
a distance to 47 Tuc of 5.0 kpc, 10% larger than previ- 
ously thought (Webbink 1985). Third, according to Howell, 
Guhathakurta & Gilliland (2000), the core radius is double 
the value used previously (De Marchi et al. 1996). In addi- 
t ion, the stellar line-of-si ght velocity dispersion at the centre 
(Meylan & Mayor 1986) is s maller than the previously ac- 



cepted value of 13.2km s"^ (Webbink 1985). The fifth up 



date is the proper motion, obtained from Hipparcos data: 
l ^a = (7.0 ± 1.0) mas yr~ ^ and ns — (— 5.3 ± 1.0)masyr~^ 



( Pdenkirchen et al. 1997] ). 

The evolutionary state of the cluster is not known pre- 
cisely. It is generally presumed (e.g., Djorgovski & King 
1984) that the core is in a pre-coUapse phase. However, other 
authors (e.g., de Marchi et al. 1996) suggest that the cluster 
underwent core collapse in the past, reaching a core radius 
similar to that of M15, and that it has now re-expanded with 
the energy provided by the creation and hardening of stel- 
lar binary systems, an im portant process in globular cluster 
evolution ( ^pitzer 1987 ). This hypothesis was based on a 
determination for the angular core radius of 12", which in- 
dicated a dynamically evolved cluster; and it can explain the 
large number of exotic objects observed in the cluste r. How- 



ever, t he more recent determina tion of core radius (Howell 



Guhatf akurta & Gilliland 2000) doubles this value, casting 
some doubt on that conclusion, since the re-expansion would 
have had to be far more vigorous. Therefore the evolution- 
ary state of the cluster remains unclear. In § ^ we discuss 
the properties of the observed binary pulsars, and their rel- 
evance for the solution of this problem is addressed in § Q 

Prior to the present work, the scarcity of coherent tim- 
ing solutions for the pulsars in 47 Tuc did not permit the 
studies that require a large sample of well-timed pulsars in 
a single cluster to be carried out. This has only been pos- 
sible to date for the eight pulsars in M15 (Anderson 1992; 
Phinney 1992, 1993). These studies of M15 provide much of 



the theoretical background for the present work and allow a 
useful comparison between M15 and 47 Tuc. 

The new observations presented here bring the number 
of pulsars with coherent timing solutions known in 47 Tuc to 
15. In § ^ we describe the observations and data reduction 
and analysis used to obtain the new timing solutions. We 
present these solutions in § ^. The new solutions lead to a 
wealth of astrophysical results which we discuss in the rest 
of the paper. In § ^ we use high-precision astrometry to 
investigate the pulsar distribution in 47 Tuc and the proper 
motion of the cluster. In § ^ we derive some constraints 
for the surface mass density of the cluster, and limits on 
the characteristic ages and magnetic fields of the pulsars. In 
§ ^ we discuss some characteristics of the binary systems, 
and present the rate of advance of periastron measured for 
47 Tuc H. Finally, in § |^, we summarise our results and 
briefly discuss the prospects for future timing observations 
of the pulsars in 47 Tuc. 



2 OBSERVATIONS AND DATA PROCESSING 

Most of the early observations of 47 Tuc were made at 
660 MHz using the 64-m radio telescope at Parkes, Australia, 
between 1989 and 1991 (Manchester et al. 1990; Manchester 
et al. 1991; for a summary of all observations see Table 
These observations led to the discovery of 10 millisecond 
pulsars. Furthe r observations of 47 T uc at 430 MHz between 
1991 and 1993 ( Robinson et al. 1995[ ) resulted in the discov- 
ery of one more pulsar, B0021— 72NQ, and phase-coherent 
timing solutions for 47 Tuc C and D, the two brightest iso- 
lated pulsars in the cluster. For the remaining nine pulsars 
then known, the paucity of detections prevented the deter- 
mination of any further coherent timing solutions. 

In August 1997, observations of 47 Tuc were resumed 
at Parkes after a four- year gap. The cluster was observed on 
126 days during the following 24 months. The majority of 
these observations were made using the central beam of a 



* Because all pulsars are detected in the same telescope beam 
pattern, and precise positions were not available originally, the 
pulsars were named B0021— 72A, B, etc. Objects A, B and K were 
later found to be spurious, so the list of pulsars in 47 Tuc starts 
with B0021-72C, now referred to as J0023-7204C, or 47 Tuc C. 
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Table 2. Summary of 10 years of observations of 47 Tuc. A^obg is the number of observing days, Vc is the central radio frequency used, 
Ai/ is the total bandwidth recorded, Afchan is the number of frequency channels per polarization, r is the most common integration time, 
tsamp is the sampling time, and tres is the time resolution for each set of observations at the central frequency. With Ai/(.han = ^^/N^^^^, 
the time resolution is obtained using t^^g = t^y^ + isampi where toM = 8-3 x (DM/cm~^pc) X {Hu^Y^^^/yUiT,) X (t'c/GHz)"^ /iS. 



Time span 


A'obs 




Ay 


-^chan 


r 


^samp 


^rcs 






(MHz) 


(MHz) 




(min) 


(ms) 


ips) 


7/89-5/92 


47 


660 


32 


128 


74 


300 


350 


5/91-11/93 


42 


430 


32 


256 


90/60 


300 


440 


1/98-4/98 


18 


660 


32 


256 


105 


125 


155 


8/97-8/99 


108 


1374 


288 


96 


280 


125 


270 



multi-beam system ( Lyne et al. 200C ) at a central frequency 
of 1374 MHz (A20 cm) with a bandwidth of 288 MHz. Incom- 
ing signals from two orthogonal polarisations were down- 
converted and filtered in a 2 x 96 x 3-MHz filter bank. After 
summing the polarisation pairs, the resulting voltages were 
one-bit sampled every 125 pis and, together with accurate 
time referencing, stored on magnetic tapes for later analy- 
sis. 

Recording the raw data in this way has two benefits. 
First, it allows off-line searches for new pulsars to be carried 
out. This strategy has been remarkably successful, essen- 
tially doubling the number of pulsars known, many of which 
are undetect able most of the tim e because of interstellar 
scintillation (Camilo et al. 2000). Second, for timing pur- 
poses, it allows us to refine the timing models iteratively, 
resulting in the ephemerides reported here. 

In the timing analysis, the raw data are de-dispersed 
and folded according to an initial ephemeris for each pul- 
sar. In the early stages, this ephemeris is determined from 
variations in the observed rotational period of the pulsar 
in the discovery and confirmation observations. It is impor- 
tant that the pulsars be detected frequently, otherwise it 
is impossible to count unambiguously the number of pulsar 
rotations between two different epochs. The 20-cm obser- 
vations made after August 1997 made use of the excellent 
sensitivity of the Parkes multi-beam system and resulted in 
a large increase in the detection rate for all the previously 
known pulsars in 47 Tuc when compared to the earlier obser- 
vations at 660 MHz. This was despite the lower flux densities 
at 1400 MHz. 

The integrated pulse profiles obtained by folding the 
data at the predicted pulse period are then cross-correlated 
with a low-noise "standard" pulse profile (see Camilo et al. 
2000 for the 20-cm profiles). This allows the determination 
of topocentric pulse times-of-arrival (TOAs) , referred to the 
observatory time standard. This was related to UTC(NIST) 
by a radio link to the Tidbinbilla Deep Space Station and 
from there to NIST by a GPS common-view system. We 
then use the tempo software package^ to calculate the 
corresponding barycentric TOAs using the assumed pul- 
s ar position and the JP L DE200 solar system ephemeris 
(tittp://ssd.jpl. nasa.gov). The differences between the mea- 
sured and predicted TOAs are used to improve the param- 
eters of the ephemeris (for further details of this process 



t http:/ /pulsar. princeton.edu/tempo 
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Table 3. Timing status for the pulsars in 47 Tuc. J2000 names 
have not been assigned to pulsars for which there is presently no 
timing solution. A horizontal line separat es the pulsars discovered 
in the early 1990s (Flobinson et al. 1995) from those reported by 
Camilo et al. (2000). A^xoA is the number of TOAs used in each 
solution, and rms is the weighted root-mean-square of the post-fit 
timing residuals in /iS. 



Pulsar 


Orbital 


Coherent timing 


solutions 




(J2000) 


solution 


Start MJD 


Final MJD 




rms 


0023-7204C 


isolated 


48494 


51406 


750 


19 


0024-7204D 


isolated 


48492 


51405 


402 


11 


0024-7205E 


known 


48464 


51406 


404 


22 


0024-7204F 


isolated 


48491 


51406 


348 


14 


0024-7204G 


isolated 


48600 


51404 


72 


26 


0024-7204H 


known 


48517 


51406 


299 


27 


0024-72041 


known 


50683 


51406 


119 


27 


0023-7203J 


known 


48491 


51406 


1114 


7 


0024-7204L 


isolated 


50686 


51406 


59 


39 


0023-7205M 


isolated 


48517 


51406 


80 


50 


0024-7204N 


isolated 


48515 


51405 


69 


14 


0024-72040 


known 


50683 


51406 


210 


20 


47 Tuc P 


known 










0024-7204Q 


known 


50689 


51405 


146 


29 


47 Tuc R 


known 










47 Tuc S 


known 










0024-7204T 


known 


50683 


51380 


125 


63 


0024-7203U 


known 


48515 


51383 


178 


20 


47 Tuc V 


not known 










47 Tuc W 


known 











see, e.g., Taylor 1992). In the early stages, the improved 
ephemeris is used to reprocess the raw data. This itera- 
tive process increases the number and quality of the TOAs, 
which in turn are used to improve the ephemeris. 

After determining the timing solutions for the 1997- 
1999 period, we re-analysed the raw data from the earlier 
observations (see Table |^. This resulted in a substantial 
increase in the number and quality of TOAs, especially for 
pulsars in binary systems, because of the much improved 
orbital ephemerides. 



3 COHERENT TIMING SOLUTIONS FOR 15 
PULSARS 

The coherent timing solutions obtained from the analysis 
described in the previous section result in a wealth of high- 
precision astrometric, spin and (for the binary pulsars) or- 
bital information. Table ^ summarises the current timing 
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Table 4. Celestial coordinates and rotational parameters for 15 pulsars in 47 Tuc at the reference epoch of MJD 51000. In this and in 
the following tables, the numbers in parentheses are l-cr confidence- level uncertainties in the last digits quoted. The positions are referred 
to the J2Q00 equinox. For all pulsars, these parameters were fitted assuming that the pulsar has the proper motion of the cluster given 
in Table |^ (see § [l:.3| ). For the DMs, the error is one or less in the last digit quoted. 



Pulsar Right Ascension Declination Period P DM 









1 m s J 




(° 


/ 


") 


(ms) 


(10-20) 


(cm pc) 


c 


00 


23 


50 


3511(2) 


—72 


04 


31 


486(1) 


5.7567799980968(5) 


—4.985(2) 


24.6 


D 


00 


24 


13 


8776(2) 


-72 


04 


43 


8323(9) 


5.3575732850382(5) 


-0.333(2) 


24.7 


E 


00 


24 


11 


1013(4) 


—72 


05 


20 


131(3) 


3.5363291476529(7) 


9.852(2) 


24.2 


F 


00 


24 


03 


8519(2) 


-72 


04 


42 


799(1) 


2.6235793491667(3) 


6.451(1) 


24.4 


G 


00 


24 


07 


956(1) 


-72 


04 


39 


683(7) 


4.040379145748(2) 


-4.215(5) 


24.4 


H 


00 


24 


06 


6989(7) 


-72 


04 


06 


789(3) 


3.210340709441(1) 


-0.162(5) 


24.4 


I 


00 


24 


07 


932(1) 


-72 


04 


39 


664(5) 


3.484992064038(1) 


-4.59(1) 


24.4 


J 


00 


23 


59 


4040(1) 


-72 


03 


58 


7720(9) 


2.1006335458586(2) 


-0.9787(4) 


24.6 


L 


00 


24 


03 


770(2) 


-72 


04 


56 


90(1) 


4.346168005784(5) 


-12.19(3) 


24.4 


M 


00 


23 


54 


485(3) 


-72 


05 


30 


72(1) 


3.676643219590(3) 


-3.832(6) 


24.4 


N 


00 


24 


09 


1835(9) 


-72 


04 


28 


875(7) 


3.053954347392(1) 


-2.186(2) 


24.6 


O 


00 


24 


04 


6492(5) 


-72 


04 


53 


751(3) 


2.6433432956679(7) 


3.032(6) 


24.4 


Q 


00 


24 


16 


488(1) 


-72 


04 


25 


149(7) 


4.033181182805(2) 


3.41(2) 


24.3 


T 


00 


24 


08 


541(2) 


-72 


04 


38 


91(2) 


7.588479792133(9) 


29.47(6) 


24.4 


U 


00 


24 


09 


8325(5) 


-72 


03 


59 


667(3) 


4.342826691451(2) 


9.524(3) 


24.3 



Table 5. Orbital parameters for the eight binary pulsars in 47 Tuc with known coherent timing solutions. For 47 Tuc H the value listed 
under Tasc is the time of passage through periastron. 



Pulsar 


Pb 


X 


^asc 


UJ 


e Lu 




(days) 


(sec) 


(MJD) 


n 




E 


2.256844818(6) 


1.981839(4) 


51000.4194595(8) 


218.5(7) 


0.000318(4) . . . 


H 


2.3576965(5) 


2.152812(7) 


51000.97359(6) 


110.494(9) 


0.070557(5) 0.059(12) 


I 


0.229792249(8) 


0.038458(9) 


51000.014865(8) 




< 0.001 


J 


0.1206649386(2) 


0.0404087(6) 


51000.0416882(3) 




< 0.0003 


O 


0.135974304(2) 


0.045157(5) 


51000.021155(2) 




< 0.0004 


Q 


1.18908405(1) 


1.462200(7) 


51000.985847(2) 


132(9) 


0.00008(1) 


T 


1.12617678(2) 


1.33847(3) 


51000.317048(3) 


55(6) 


0.00037(4) 


U 


0.4291056829(6) 


0.526953(5) 


51000.0705011(6) 


341(7) 


0.00015(2) 



status for the pulsars known in 47 Tuc. The orbital param- 
eters for 47 Tuc S and T were first determined using a new 
technique described by Freire, Kramer & Lyne (2001), while 
the orbital parameters of the remaining binaries had been 
determined earlier (Robinson et al. 1995; Camilo et al. 2000). 
The timing solutions were obtained as follows. 

Because the dispersion in proper mot ions among pul- 
sars is expected to be very small (see § 4.3) we assume that 



all pulsars have the average proper motion of the cluster, 
whose value we take to be that determined from Hipparcos 
data (Table [l|). Using tempo, the available 1991-1999 TO As 
were fitted to a model for each pulsar containing celestial co- 
ordinates, spin parameters, and binary elements where rel- 
evant. Dispersion measures (DM) were obtained separately 
by measuring frequency-dependent delays across the 288- 
MHz bandwidth available at 1400 MHz, and the resulting 
values coincide with those given by Camilo et al. (2000) . We 
also fitted astronomically meaningless time offsets between 
groups of TO As obtained at different frequencies and with 



different time resolution. We do this because it is difficult to 
make a proper absolute alignment of pulse profiles obtained 
at different frequencies owing to variations in pulse shapes. 

In Table ^ we present the positions, rotational param- 
eters, and DMs obtained in this way for 15 pulsars. The 
corresponding timing residuals for the seven isolated pul- 
sars are displayed in Fig. |^ as a function of time. 

The orbital elements for the eight binary pulsars with 
timing solutions are presented in Table ^ while the corre- 
sponding timing residuals are shown in Fig. ^ as a function 
of time, and in Fig. ^ as a function of orbital phase. 

For each pulsar we list the five Keplerian parameters: bi- 
nary period Pb, projected semi-major axis fight travel time 
X, time of passage through the ascending node Tasc and, 
where measurable, longitude of periastron uj and eccentric- 
ity e. For 47 Tuc H we indicate the time of passage through 
periastron, rather than through ascending node, and we also 
list the measured rate of advance of periastron lj. For low- 
eccentricity binaries (all but 47 Tuc H), in which some or- 
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Figure 1. Post-fit timing residuals as a function of date for seven isolated pulsars. In this and in Figs. Q and |^ the crosses, circles 
and filled diamonds represent residuals at 1400, 660 and 430 MHz respectively. At right we display corresponding histograms of timing 
residual/uncertainty for individual TOAs. 
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1991 1992 1993 1994 1998 1999 -2 2 



Date Residual/a. 

Figure 2. Post-fit timing residuals as a function of date for eight binary pulsars and corresponding histograms of timing resid- 
ual/uncertainty for individual TOAs. 
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0.5 1 



Orbital phase 

Figure 3. Post-fit timing residuals as a function of orbital phase for the eight binary pulsars whose residuals are shown as a function 
of date in Fig. ^. Phase is measured relative to periastron for 47 Tuc H and to ascending node for the remaining pulsars. 47 Tuc O is 
eclipsed for a portion of its orbit, and we do not include TOAs from this region in the timing fit. 
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Table 6. Proper motions of five pulsars in 47 Tuc. The average is 
obtained by weighting the proper motions by the inverse of each 
uncertainty. 



Pulsar 


Mo, 

(mas yr"-*-) 


A»4 

(mas yr~^) 


C 


7.8(1.6) 


-3.2(1.0) 


D 


7.7(1.8) 


-2.6(1.6) 


E 


9.1(2.3) 


-3.9(2.1) 


F 


6.6(2.6) 


-3.7(2.0) 


J 


4.6(0.8) 


-3.7(0.8) 


Average: 


6.6(1.9) 


-3.4(0.6) 



bital parameters are highly covariant in standard binary fits, 
we used the binary model ELLl as implemented in tempo 
to determine the solutions presented (see Lange et al. 2001). 

Finally we attempted to measure proper motions for 
those pulsars where we had a combination of a long time 
baseline and generally high-quality data. We did this with a 
straightforward tempo fit, and list the results in Table ^ for 
five pulsars. We also determined proper motions by compar- 
ing the positions of the pulsars in 1992-1993 and 1998-1999: 
for each of the two independent position measurements, no 
proper motion was assumed, and all TOAs used at each 
epoch had the same frequency (430 MHz in early 1990s, 
1400MHz in late 1990s), thus avoiding any problems with 
the alignment of standard profiles at different frequencies. In 
these fits we used the accurate rotational and binary param- 
eters obtained previously, so only the celestial coordinates 
were determined. A proper motion was then calculated from 
the differences in position between the two epochs. These 
measurements are consistent with those made in the global 
fits, but they are less precise. 

All uncertainties presented in Tables ^|-^ are our best 
estimates of realistic 1-cr confidence levels, and in most cases 
are twice the formal fit uncertainties from tempo. 

The main factors limiting the precision of the timing so- 
lutions presented in this paper are the relatively low signal- 
to-noise ratio of the pulsars and the low detection rate, par- 
ticularly in the observations of the early 1990s. For some of 
the pulsars, the 20-cm fiux densities are so low that they are 
detectable only in about 10% of all observations (see Camilo 
et al. 2000). Additionally, the data obtained between July 
1989 and May 1991 were not used in any of the fits be- 
cause the difference between the Parkes clock and UTC is 
not known to better than 50 /is. This introduces significant 
errors that would compromise several of the measurements. 



Table 7. East-west (x) and north-south (y) offsets of 15 pulsars 
from the centre of the cluster, assumed here to be exactly at 
a = 00''24™05!29 and S = -72°04'52'.'3. The errors in the off'sets 
are 2 or less in the last digits quoted. The angular distance of each 
pulsar from the centre of the cluster, 8j_, has an accuracy limited 
by the uncertain absolute position of the cluster (Table 111). 



Pulsar 


X 


y 


0± 




(arcmin) 


(arcmin) 


(arcmin) 




1 1 ^41 
— i . 10^1 




1 91 


D 


0.6634 


0.1411 


0.68 


E 


0.4489 


-0.4639 


0.65 


F 


-0.1111 


0.1583 


0.19 


G 


0.2060 


0.2103 


0.29 


H 


0.1088 


0.7585 


0.77 


I 


0.2041 


0.2106 


0.29 


J 


-0.4547 


0.8921 


1.00 


L 


-0.1174 


-0.0767 


0.14 


M 


-0.8347 


-0.6404 


1.05 


N 


0.3008 


0.3904 


0.49 





-0.0495 


-0.0242 


0.06 


Q 


0.8651 


0.4525 


0.98 


T 


0.2512 


0.2232 


0.34 


U 


0.3509 


0.8772 


0.94 




X (arcmin) 



4 PULSAR POSITIONS 

The timing solutions presented in the previous section bring 
the number of pulsars in 47 Tuc with accurately known po- 
sitions to 15, and with measured proper motions to five. For 
easy reference, we derive from the coordinates in Table ^ 
the east-west (x) and north-south (y) offsets of the pulsars 
relative to the cluster centre. These are presented in Table ^ 
and depicted in Fig. ^ In the remainder of this section, we 
discuss the implications of these positional measurements. 



Figure 4. Positions of the 15 pulsars in 47 Tuc with a timing 
solution plotted in the piano of the sky. The positions are given 
as east-west (x) and north-south (y) offsets from the centre of the 
cluster (see Table |7 ) . The pulsars indicated with © have positive 
observed period derivatives (P), while those indicated with 
have negative P. The central circle indicates the core radius. 
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9^ (core radii) 

Figure 5. Cumulative radial distribution of the pulsars in 47 Tuc. 
This is compared with the expectation of a distribution of the type 
n(r) oc with a cutoff at 4 core radii (solid line) and 3 and 5 
core radii (dotted lines). 




Figure 6. Globular cluster potential as a function of distance 
from the centre given by a King model of 47 Tuc. The maximum 
squared pulsar velocity at the centre is indicated by the height of 
the dashed line above W{r = 0) = — 8.7d^(0). The pulsars ca nno t 
travel further than about 5 core radii from the centre (see § 



4.1 Radial distribution 

The radial distribution for the pulsars in 47 Tuc was first 
addressed by Rasio (2000), based o n an earlier versi on of 
some of the results we present here (Freire et al. 2000). We 



now outline the main characteristics of this distribution and 
draw some conclusions from it. 

The most striking characteristic of the pulsar distribu- 
tion is that all the pulsars in 47 Tuc with measured positions 
lie within ll2 (3 core radii) of the centre of the cluster, de- 
spite the fact that the tidal radius of 47 Tuc is about 40' 
(Table |i|). This distribution, as we shall now see, is not an 
artifact introduced by the size and shape of the Parkes tele- 
scope radio beam pattern, which has a half-power radius of 
~ 7' at 1400 MHz. 

The flux densities for most of the pulsars have been 
calculated by Camilo et al. (2000) . Among the pulsars with 
a known solution, 47 Tuc N has the lowest fiux density (0.03 
± 0.01 mjy). Supposing that this value is the lower limit for 
the flux density of a pulsar for which we can obtain a timing 
solution, weak pulsars like 47 Tuc U, with a flux density of 
0.06 ± 0.01 mJy, are detectable in a circle with a radius of at 
least 7'. There are 10 pulsars with at least this flux density, 
and the fact that none of these is seen outside a radius of 
l!2 is therefore a true feature of the pulsar distribution, and 
not an artifact due to the shape of the beam. It should also 
be noted that the original surveys at 660 MHz (Manchester 



et al. 1991), with a larger telescope beam, covered an area 



approximately 100 times the size of the roughly arcmin-scale 
central region in which we now localize the 15 pulsars. 

The present pulsar distribution is essentially an equilib- 
rium distribution. The characteristic ages of these pulsars 
are of the order of 10^ yr, while the relaxation time (the 
time it takes for a pulsar to change kinetic energy signifi- 
cantly through interactions with other stars) in the core of 



47 Tuc is of the order of 10** yr ( pjorgovski 1993[ ). There- 
fore, the pulsars are expected to be in thermal equilibrium 
with the remaining stars, and their distribution should be 
merely a function of their mass and of the potential of the 
cluster. Indeed, the de-projected pulsar distribution is con- 
sistent with n(r) ~ until about 4 core radii, decreasing 
sharply outside this radius (see Fig. 

This concentration of pulsars near the centre of the clus- 
ter is a consequence of mass segregation: if all objects in a 
cluster have reached thermal equilibrium, i.e., all stellar pop- 
ulations have a similar average kinetic energy, then the most 
massive populations, among which are the pulsars, will have 
smaller velocities, and therefore will dwell deeper in the po- 
tential well of the cluster. To quantify this, we note that the 
kinetic energy of the lighter stars in the cluster c ore, with 
mass nims , is limited by the escape velocity Ve. = ^y—2W(0). 
Equating the resulting limit of kinetic energy to that of the 
pulsars results in a limit to the velocity of pulsars bound in 
the cluster core; 



(1) 



where nip is the pulsar mass. As pulsars move away from 
the centre of the cluster under the influence of the cluster's 
gravitational potential their velocities decrease, eventually 
reaching zero at a maximum possible distance rum from the 
centre of the cluster. 

The potential energy U{r) at is calculated from 
conservation of energy: 

[/(riim) = mpW{riira) = mpW{0) + (2) 

where W(r) is the cluster potential. Using equation |l] we 
obtain 



W^(nim) = M^(0) 1 



(3) 
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This limit, with rrims = 0.8 Mq and rrip = 1.45 M© (the 
average mass for isolated and binary pulsars) , is represented 
by a dashed line in Fig. ^ where we use a King model for 
the gravitational potential of 47 Tuc (King 1966; Phinney, 
priv. comm.). From this we derive rum ~ 5rc, which is very 
close to the inferred 4 rc. Therefore, the confinement of the 
pulsars in the inner region can be explained, to first order, 
as resulting from the shape of the potential and the upper 
limit for the kinetic energy of the pulsars. The fact that all 
pulsars are located close to the core is a constraint on any 
future more accurate model of the gravitational potential of 
the cluster. 

Spitzer (1987) demonstrates that the distributions of 
two stellar species of mass rrii and rrij in thermal equilibrium 
with all stars in a cluster are related by 

ni{r) (X nj{r)™^^"^^ . (4) 

Therefore, two stellar populations with similar masses in 
thermal equilibrium should have similar radial distribu- 
tions. Some of the most massive blue stragglers know n in 
47 Tuc have masses of 1.3-1.6 Mq ( iGilhland et al. 1998|). Be- 
cause t his is the mass range expected for pulsars (Thorsett 
& Cha] :rabarty 1999), incl uding binaries such as those in 
47 Tuc (Camilo et al. 200C), both populations should have 
similar radial distributions, which is apparently the case 



( Rasio 200C ) . However we caution that the population con- 
sidered by Rasio (2000) is not altogether homogeneous, hav- 
ing a relatively wide range of masses. Ultimately it will be 
useful to compare the radial distribution of just the most 
massive blue stragglers with that of pulsars, and in partic- 
ular to determine whether they display the abrupt radial 
cutoff observed for pulsars. 

Lastly, Spitzer shows that outside 3 core radii, the mass 
distribution for the dominant stellar species should be of the 
type n(r) oc . Knowing this and the fact that the radial 
distribution of the pulsars in M15 is n{r) oc r~^'^ (Anderson 
1992 |), ij 'hinney (1992) used equation^ to infer that the dom- 
inant stellar species near the centre of M15 is 0.9 Mq white 
dwarfs. In principle we could use the same argument to show, 
based on the observed pulsar distribution in 47 Tuc, that the 
dominant stellar species in the central region of this cluster 
is made of 1.3-1.6 Mq objects. However, the fact that the 
pulsars lie within (and not outside) 3 core radii may compli- 
cate this interpretation (although the "pulsar core" should 
be more compact than that of the cluster). 

4.2 A triple system with two pulsars? 

There is a remarkable coincidence between the celestial co- 
ordinates of 47 Tuc G, an isolated pulsar, and 47 Tuc I, a 
binary pulsar (see Table The projected angular distance 
between these two pulsars is only 120 mas. 

To estimate the chance probability of finding two pul- 
sars in such close proximity, suppose that pulsars are ran- 
domly distributed in a disk of radius 18" (the projected an- 
gular distance of these two pulsars from the centre of the 
cluster) . For a given pulsar, the probability of finding a sec- 
ond one at a projected distance of 0'.'12 or less is proportional 
to the area of a disk of radius 0'.'12 divided by the area of a 
disk of radius 18" or ~ 5 x 10~^. The probability of being 
at a distance larger than 0"12 is 0.99995. For the five pul- 
sars known within 18", we have 10 difi'erent possible pairs. 
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Figure 7. Measured proper motions of five pulsars in 47 Tuc. The 
semi-axes of the ellipses represent the 1-cr errors for the proper 
motions quoted in Table ^. The v alue for the proper motion of the 
cluster (Odonkirchen et al. 1997) is indicated by the thick circle. 



The filled ellipse indicates the motion of the cluster and its error 
deduced by averaging the measured pulsar proper motions. 



The probability of not finding any pulsar within 0'.'l2 of any 
other pulsar is 0.99995^° = 0.9995, i.e., the overall probabil- 
ity of finding one or more pulsars projected within 0'.'12 of 
any other is 5 x 10~*, about one in 2000. 

If the pulsars are near to each other but not interact- 
ing significantly, then the acceleration caused by the cluster 
should be similar for both of them. In § ^ we see that their 
values of P/P are similar, with a difference of 3 x 10~^* s~^, 
in a possible range ~ 20 times larger. Therefore, the com- 
bined probability of finding the pulsars with such close 
projected separation and with such similar accelerations is 
about one in 40,000. 

Given this low formal probability of chance coincidence, 
it is worth considering that these two pulsars may be in 
a hierarchical triple system with a major axis of at least 
600 a.u. A system with two components 600 a.u. apart is not 
likely to remain bound for long in the dense environment of 
a globular cluster. By considering the stellar flux with one- 
dimensional speeds of Vz{r), an estimate for the mean time 
a binary system can survive in a dense environment (i.e., 
before it gets hit by another cluster member) is 



(5) 



In this expression n(r) is the density of stars, Vz{r) is the 
one-dimensional stellar velocity dispersion near the binary 
and a is the cross-sectional area for the interaction of the 
binary with a star from the cluster. The King model for 



47 Tuc predicts n(r 
-1 



4 X 10" pc"^ for this location. With 
Vz{r) ~ lOkms"^ we estimate r ~ 10''-10^ yr. This time is 
comparable to the orbital period of the putative binary sys- 
tem, and 10~^-10~^ the age of the cluster. Thus, despite the 
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low formal probability of a chance coincidence, this simple 
calculation argues against the reality of such an association. 
Ultimately, the reality or otherwise of this system can be 
tested by future measurements of higher period-derivatives 
for the two pulsars. 



4.3 Proper motions 

The proper motions listed in Table |^ are displayed in Fig. ^ 
together with the value measured for the proper motion of 
the cluster by Odenkirchen et al. (1997). At the 2-a level all 
these proper motions are consistent with each other. 

If a pulsar is moving with a transverse velocity Vp 
(kms~^) relative to the centre of the cluster, the difference in 
proper motions is A/x = 0.21Wp/_Dmasyr~^, where D is the 
distance to the cluster in kpc. The escape velocity of 47 Tuc 
is about 58 kms~^ (Table hi). Therefore, in order to detect a 
pulsar in an escape trajectory at the 3-a level, we would have 
to measure relative proper motions of about 2.4masyr~^ 
with a precision better than 0.8masyr~^, and we are ap- 
proaching this level for 47 Tuc J (Table ^). If the pulsars are 
in thermal equilibrium with the surrounding stars, their rel- 
ative motions are of the order of the central velocity disper- 
sion or less, i.e., ~ 12kms~^. Thus, detecting their relative 
motions implies measuring relative proper motions of about 
0.5masyr~^ with a precision of 0.15masyr~^. The typical 
precision for the positions of the weakest pulsars is around 
5mas. In order to obtain a precision of 0.15masyr~^ for 
these we would have to measure the positions of the pulsars 
again in ~ 50 years' time. However, for the brightest pulsars, 
with positional accuracies of about 1 mas, it will suffice to 
make another measurement in about 10 years, using current 
data-acquisition systems. 

The assumption of thermal equilibrium for the pulsars 
implies that, with current timing precision, a measurement 
of the proper motions of these pulsars is in effect a measure- 
ment of the proper motion of the cluster. We can therefore 
make a weighted average of the five proper motions and 
determine the motion of the cluster, knowing that the pe- 
culiar motions of the pulsars should be smaller than the 
errors in the individual measurements. The weights chosen 
are simply the inverse of the uncertainties in each coordi- 
nate. For the currently-observed proper motions, we find the 
weighted average value to be fia = (6.6 ± 1.9) masyr"^ and 
fj,s = (—3.4 ± 0.6) masyr"^. The uncertainties are taken to 
be the sum in quadrature of the weighted dispersion of the 
values of proper motion for the pulsars about the average, 
and 0.5masyr~^, which accounts for the expected actual 
dispersion of proper motions. 



5 PULSAR ACCELERATIONS IN THE 
CLUSTER POTENTIAL 

Nine of the 15 pulsars in Table ^ have negative period deriva- 
tives (P). This has been observed before for some pulsars 
located in globular clusters (e.g., Wolszczan et al. 1989). 
Rather than being due to intrinsic spin-up, negative period 
derivatives are thought to be caused by the acceleration of 
the pulsar towards the Earth in the cluster potential (see 
Fig.l). 

In § 5.1 we use the observed period derivatives of the 




. ■ \/' -. j" • To Earth . ;, . 

Figure 8. Geometrical parameters used in § 11 is a plane per- 
pendicular to the line of sight that passes through the centre of 
the cluster, at point O, with core radius r^. For a pulsar at point 
P, R± is the projected distance to, and a is the acceleration of the 
pulsar towards, the centre of the cluster. The line-of-sight compo- 
nent of a, ai, is the only one detectable from the Earth (located 
in the plane of the figure). 



pulsars to estimate lower limits for the surface mass density 
of the cluster. We then consider the aptness of a King model 
together with accepted parameters for 47 Tuc in describing 
the accelerations experienced by the pulsars as inferred from 
the period derivatives (§ 5.2). We derive limits for the ages 
and magnetic fields of the pulsars in § 5.3. 

The observed period derivative, Pobs, is the sum of the 
pulsar's intrinsic spin-down. Pint, and the effect of the accel- 



eration along the line of sight, a;. This sum may be negative 
if a negative a; contribution is not exceeded by a positive 
Pint; i.e., a negative Pobs implies \ai/c\ > P/Pnt- It also 
implies \ai/c\ > |P/Pobs|. 

The acceleration along the line of sight has several com- 
ponents, denoted by as, ac, and ac in 



as ac ac 
P / c c c V P 

obs 



where the terms represent: 



(6) 



1. The centrifugal acceleration (Shklovskii 1970), given 



by 

as_ 

c 



m'p> 



(7) 



This term, using the values in Table |^, amounts to as/c = 
(+9 ±3) X lO^^^s^V 

2. The difference in Galactic acceleration along the line 
of sight between a given object and the bary centre of the 
solar system (ac) - This is a function of the object's Galactic 
latitude b, longitude /, and distance from the solar system 
D. For 47 Tuc we have b = -44.9° , / = 305.9° and D = 5.0± 
0.4 kpc. Using Paczynski's (1990) model of the gravitational 
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Table 8. Constraints on projected surface mass density for 4 7 Tu c 
obtained from nine pulsars with negative (-P/P)obs (see § 5.1). 
The distance to the cluster is taken to be 5.0 kpc. 



Pulsar 


{P/P)ohs 

(10-1'' s-1) 


(pc) 


(10=* Mq pc-2) 


M'cyi(< 6»x) 
(103 M©) 


C 


-0.91 


1.75 


> 18 


> 170 


D 


-0.11 


0.99 


> 2.2 


> 6.6 


G 


-1.09 


0.43 


> 21 


> 12 


H 


-0.10 


1.11 


> 1.9 


> 7.5 


I 


-1.36 


0.43 


> 27 


> 15 


J 


-0.51 


1.46 


> 10 


> 67 


L 


-2.85 


0.20 


> 56 


> 7.3 


M 


-1.09 


1.53 


> 21 


> 160 


N 


-0.76 


0.72 


> 15 


> 24 



potential of the Galaxy, we obtain aa/c = (—4.5 ± 0.2) x 
10-1^ s-i. 

3. Accelerations due to the gravitational field of the 
globular cluster and its individual stars (ac)- This is the 
most interesting contribution from an astrophysical point of 
view. 

Contributions 1 and 2 total (+5±3) x 10"^^ s~^. Hence- 
forth, (P/P)obs indicates the measured value of P /P minus 
these contributions, and a; refers solely to ac- AH conclu- 
sions in this section apply only to the pulsars with known 
timing solutions. 



5.1 Projected surface mass density of 47 Tuc 

In Fig. ^ are shown some useful geometrical parameters we 
use in the remainder of this section. The plane 11 passes 
through the centre of the cluster and is perpendicular to 
the line of sight, and the core radius is given by rc = D6c- 
A particular observed pulsar has unknown distance to the 
centre of the cluster, r, and to 11, Z — we know only the 
projected distance of the pulsar to the centre of the cluster, 
R±_ — D9±. The acceleration along the line of sight, ai, is 
the only component potentially detectable from the Earth. 

We now derive constraints for the surface mass density 
of the cluster in its inner regions. According to Phinney 
(1993), the following expression is valid to within ~ 10% for 
all plausible cluster models, and for all 6±: 



5.1 X 10" 



1.1 



GMcyl(< 6>i 

E(< Si) 



(l03 Me pc- 



(8) 



where Mcyi(< 0_l) is the mass of all the matter with a pro- 
jected distance smaller than 6±, and E(< 6'x) is the cor- 
responding projected surface mass density. The limit on S 
does not depend on estimates of the distance to the clus- 
ter or core radius. The pulsars with negative observed pe- 
riod derivatives provide a lower bound on ainiax(^x) and 
therefore, with equation ^, on S and Mcyi- These limits are 
presented in Table ^. 



5.2 Accounting for the pulsar accelerations 

We have used a King model of 47 Tuc to calculate the grav- 
itational potential as a function of distance to the centre of 
the cluster, W{r/rc), divided by Uz(0)^ (Fig. |). The input 
for such a model is a single parameter, the logarith m of the 
core radius divided by the tidal radius ( King 196^ . 

The radial derivative of W{r /rc)/vz{QY yields a "nor- 
malized acceleration", an{r-/rc)- In order to calculate accel- 
erations j4 in m s~'^, we need two more cluster parameters: 
the central dispersion of line-of-sight velocities, Wz(0), and 
the core radius, rc = DOc- A is obtained from 



A{r/rc) = 



anir/rc)vi{0) 
DOc 



(9) 



To calculate the contribution to (P/P)obs of the accel- 
eration of a pulsar at a given radius we multiply A{r /rc) by 
l/r, so as to obtain the line-of-sight component of the accel- 
eration. We do not, of course, know the actual distance r of 
each pulsar from the cluster centre, only its projection R±. 
We will therefore compare the "observed acceleration" of 
each pulsar to the ensemble average along the relevant lines 
of sight calculated from the King model, taking the pulsar 
distribution in the cluster to have a density proportional to 

and a cut-off radius of 4 core radii (see § [4.l| ). 

The average acceleration along each line of sight (inte- 
grated over only one hemisphere) is thus given by 



Aa{Rl 



1 



— A{r/rc)-dl, 



(10) 



where A is a normalizing constant such that A~^/r'^ is the 
local linear density of simulated pulsars, and N is their total 
number in a column along the integration path. 

How do these "average" theoretical accelerations along 
each line of sight compare with the observed accelerations? 
The observed acceleration, Ao, can be obtained for each pul- 
sar by subtracting the intrinsic P/P from {P/P)obs- Un- 
fortunately, the intrinsic period deriva tive is not known. 

where we also 



5.3 



We derive limits on this quantity in 
find that the average characteristic age of the pulsars is 
greater than ~ 1 Gyr. Here we assume that each pulsar has 
a characteristic age of 3 Gyr and subtract the corresponding 
P/P — +0.5 X lO"^'^ s"'^ from each observed value (listed in 
Table |l^) to obtain an estimate for Ao- These, along with 
values of Aa calculated from equation ^ are presented in 
Table |. 

We are now in a position to compare the "observed" and 
average predicted accelerations. We do this for each pulsar 
in the last column in Table ^ This value is significantly 
biased for some individual pulsars by the subtraction of a 
fixed P/P — 0.5 X 10~^^s~^, our first-order attempt at ac- 
counting for the intrinsic period derivatives. However, the 
average of all values in the last column of the table should 
be approximately unity, if the modeling described above is 
to be consistent with the actual accelerations experienced 
by the pulsars — which indeed it is, despite considerable 
uncertainty: \Ao/Aa{R±)\ = 1.2 ± 0.7. 

As an additional consistency check on the King model of 
the cluster potential we performed a Monte Carlo simulation 
of the pulsar population of 47 Tuc. The aim of this simula- 
tion is to see, given reasonable input assumptions, whether 
we can reproduce the distribution of (P/P)obs. 
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Table 9. Comparison of estimates of "observed" pulsar accelerations Ao with predicted aver age accelerations. Aa{R±), derived using a 
King model for 47 Tuc with D = 5.0±0.4kpc, Vz{0) = 11. 6± 1.4kms~^ and 9c = 23'.'1 (see § |5.2| ). The maximum acceleration predicted 
by the model for each projected pulsar offset from the centre of the cluster, used in Fig. 11 is also listed. 



Pulsar 


\Ao/c\ 
(lO-i'^s-i) 


\Aa{R_L)/c\ 

(lO-i'^s-i) 


Wl max(-R_L)/c| 

(10-1^ s-1) 


\Ao/Aa{R_L)\ 


C 


1.38 


0.6 ±0.2 


1.0 ±0.3 


2.4 ± 1.0 


D 


0.57 


1.3 it 0.5 


1.8 ± 0.7 


0.4 ± 0.2 


E 


2.27 


1.3 it 0.5 


1.9 ± 0.7 


1.7 ± 0.7 


p 


1 


9 1 -1- n s 
z. 1 in u.o 


o.u m i.o 


n Q -t- n A 




1 56 


2 it 0.7 


3.2 ± 1.1 


0.8 ± 0.3 


H 


0.56 


1.1 it 0.4 


1.6 ± 0.6 


0.5 ±0.2 


I 


1.83 


2.0 it 0.7 


3.2 ± 1.1 


0.9 ±0.4 


J 


0.98 


0.8 it 0.3 


1.2 ± 0.5 


1.2 ±0.5 


L 


3.32 


2.0 it 0.7 


3.8 ± 1.4 


1.6 ±0.7 


M 


1.55 


0.8 it 0.3 


1.2 ± 0.4 


2.0 ±0.9 


N 


1.23 


1.6 it 0.6 


2.4 ±0.8 


0.8 ±0.3 


O 


0.63 


1.4 it 0.5 


4.0 ± 1.4 


0.5 ±0.2 


Q 


0.33 


0.9 it 0.3 


1.3 ± 0.4 


0.4 ±0.2 


T 


3.37 


2.0 it 0.7 


3.0 ± 1.0 


1.7±0.7 


U 


1.68 


0.9 it 0.3 


1.3 ± 0.5 


1.9 ±0.7 




Figure 9. Cumulative distributions of (P / P)ohs./Wl max/c| for 15 
observed pulsars, and for a simulated populatio n (s mooth curve), 
normalized by the observed distribution (see § |5.2|). 



In the simulation, we generate a spherically symmet- 
ric population of 10^ pulsars with a radial distribution of 
the type n(r) oc and a cutoff of 4 core radii. Each 
of the pulsars was randomly assigned an age from a flat 
distribution ranging between 500 Myr and 10 Gyr — the 
mean P/P = 0.5 x 10~^'^s~^ corresponding to what was 
assumed before. These assigned ages correspond to the in- 
trinsic P/P of each of the model pulsars. Finally, given the 
position of each of the model pulsars relative to the clus- 
ter centre, we calculate the contribution to P/P from the 



acceleration in the cluster determined from the King model 
with D = 5.0 kpc, w^O) = 11.6 kms"^ and Oc = 23'.'l. By 
adding this contribution to the intrinsic P/P we have, for 
each pulsar, a model observed P/P which can then be di- 
rectly compared to the sample of 15 pulsars for which we 
have timing solutions. 

The results of this simulation are shown in Fig. ^ where 
we present a cumulative plot of {P / P) ohs / \o.i mayL{R±) /c\ for 
the real and simulated samples. Given the relatively small 
size of the observed distribution, and the straightforward 
simulation that we have performed, the agreement between 
the model and observed samples is good. 

Our overall conclusion is therefore that a King model 
for 47 Tuc, with D = 5.0±0.4kpc, v^O) = 11.6±1.4kms-\ 
dc = 23'.' 1 and a small contribution from the intrinsic P/P 
of the pulsars, provides a good description for the observed 
values of P/P. 



5.3 Limits on pulsar ages and magnetic fields 

We now obtain one-sided limits for the characteristic age 
and the inferred surface dipole magnetic field strength of the 
pulsars, by calculating the maximum intrinsic period deriva- 
tive for each pulsar. Pint max. First, suppose we have a reli- 
able estimate for the maximum possible acceleration for each 
pulsar's line of sight, a;niax(-R±)- Next, assume that each 
pulsar actually experiences such a maximum (negative) ac- 
celeration. Under such conditions, an observed \P/P\ that is 
different from |a;max(P±)/c| is caused by a finite (positive) 
intrinsic P/P. Hence, (P/P)int < |a,n.ax(P±)/c| + (P/P)obs. 

After determining Pint max we derive a minimum char- 
acteristic age of the pulsar, Tc > P/(2Pintmax), and 
an upper limit for the magnetic field, B < 3.2 x 

10l^(PPint max) '/'Gauss. 

Only a reliable upper limit for the acceleration along 
each line of sight remains to be calculated. Phinney (1993) 
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Figure 10. The values for (P/P)obs (Tabic 0) are plotted ver- 
sus d x_ for the 15 pulsars with coherent timing solutions (stars). 
We also plot the maximum acceleration along each line of sight 
calculated from a King model with the nominal distance and cen- 
tral velocity dispersion (solid line), and accounting for the uncer- 
tainties in those quantities (dashed lines). For all the curves and 
pulsarpositions, we assumed dc to have its nominal value. See 
Table 1 for cluster parameters and their uncertainties. 



straints for a; max(i?±), more stringent limits can be derived. 
For this reason, we now compare the observed pulsar P/P 
with the maximum accelerations along each line of sight cal- 
culated from the King model (Fig. [10|). The limits for the 
pulsar parameters derived from this model are also presented 
in Table |lo[ As expected, these are more constraining than 
those derived before. In particular, the characteristic ages 
for 47 Tuc C and M exceed 3 Gyr, and a simple average of 
the lower limits on age for all pulsars yields Tc > 0.9 Gyr, 
while the magnetic field strengths of 47 Tuc J and M are 
less than 3 x 10* Gauss, very close to the lowest values ob- 
served among Galactic disk pulsars. These values depend, of 
course, on the particular cluster model used in calculating 
the accelerations. 

Limits for the characteristic ages and magnetic fields of 
the pulsars in M15 have been derived by Anderson (1992), 
using a mass model for the cluster based on the surface lumi- 
nosity density and assuming a constant mass-to- light ratio. 
The majority of pulsars in th at cluster have characteristic 
ages ~ 10^" yr (Anderson 1992), but there are signs of a rela- 
tively recent burst in pulsar formation. At least two pulsars 
have maximum characteristic ages ~ 10* yr and magnetic 
fields of about (10-20) x 10^ Gauss. This is attributed to the 
ongoing core collapse in M15, leading to much increased re- 
cycling of pulsars. In 47 Tuc, on the other hand, there are 
no visible signs of such a recent burst in pulsar formation: 
all the known pulsars for which a timing solution has been 
obtained have characteristic ages (probably much) greater 
than 170 Myr, weak magnetic fields, and very short rota- 
tional periods. 



calculated the maximum acceleration expected near the cen- 
tre of a globular cluster. 



Oil 



- o 



(11) 



This expression is valid to within ~ 10% for R± ^ 2rc, 
provided that the function Vz{R±) is accurately known. 

The values of P/P observed for the pulsars (reflecting 
in part the gravitational potential of the cluster) do not 
appear to decrease greatly with R±_ (Fig. so we first as- 
sume conservatively that aimax(Px) does not decrease with 
R±. Therefore, the absolute upper limit for aijnax(P±) is 
Hi max (0). Within the constraints imposed by the parame- 
ters and uncertainties in Table |l], the maximum possible 
value for a; max at the centre of the cluster is obtained with 
D = 4.6 kpc, v^O) = 13.0 kms"^ and 9c = 2l'.'4: according 
to equation |ll|, ai„^B,^/c = 5.74 x lO"'^^ s~^. 

The limits for the pulsar parameters obtained with this 
constant maximum acceleration are presented in Table |lo[ 
All characteristic ages are larger than 170 Myr and all pul- 
sars have magnetic fields lower than 2.4 x 10^ Gauss. These 
values are consistent with those measured for Galactic mil- 
lisecond pulsars (e.g., Camilo, Thorsett & Kulkarni 1994). 

Note that these limits are extremely conservative, and 
are independent of any detailed modeling of the cluster. It 
is nevertheless clear that even with such a crudely over- 
estimated aimax(P±) we can derive useful limits for the pa- 
rameters of the pulsars. 

By modeling the cluster and obtaining better con- 



6 BINARY PULSAR SYSTEMS IN 47 TUC 

The binary systems containing pulsars in 47 Tuc are di- 
vided into two main g roups, segregated by companion mass 
( [Camilo et al. 2000| ). The first is composed of binaries 
with very short orbital periods (1.5-5.5 hr) and compan- 
ion masses ~ 0.03 M©. We have timing solutions for three 
of these (47 Tuc I, J and O). The second group has or- 
bital periods in the range 0.4-2.3 d and companion masses 
~ 0.2 Mq; we refer to this group as "normal" binaries. Of 
these, 47 Tuc E, H, Q, T and U have a timing solution at 
the moment. 

Coherent timing solutions give extremely detailed or- 
bital information, and often allow the measurement of very 
small eccentricities. Most millisecond pulsar-white dwarf 
systems in the Galactic disk have very low, but measur- 
ably significant eccentricities (e.g., Camilo 1999). These are 
thought to result from tidal interactions between the neu- 
tron star and convective cells in the envelope of its compan- 
ion during the giant phase of its evolution, and a relation- 
ship between orbit al period and eccentricity can be derived 
for such systems (Phinney 1992). The millisecond pulsar- 



white dwarf systems in globular clusters often ha ve larger 
eccentricities than those in the disk of the Galaxy ( Phinney 



1992). These are thought to result from interactions with 



other stars in the cluster: the denser the environment, the 
more likely one is to find a relatively large degree of eccen- 
tricity in the system. 

The five eccentricities measured for binary pulsars in 
47 Tuc (for the normal binaries 47 Tuc E, H, Q, T and U; 
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Table 10. Limits on period derivative, characteristic age and magnetic field for the pulsars in 47 Tuc. The {P/P)ohs have been corrected 
for the acceleration terms as and ac using D = 4.6 kpc. The values annotated with "2" were derived from the maximum accelerations 
(also indicated) calculated with an isotropic single-mas s Ki ng model for 47 Tuc with the same value for the cluster distance, Vz{0) = 
13.0kms~^ and the nominal value for core radius (see § 5.3). The values annotated with "1" were calculated with the same parameters, 
except for angular core radius, where the lower limit, 6c = 21'.'4, was used. In this case li max (Kx ) /c = (iimax(0)/c = 5.74 X lO^^'^s"^ 
(eq. §). 



Pulsar P (P/P)ohs Antl Tcl Bl |a;n,ax/c| -Pint 2 Tc 2 B2 

(ms) (10-17 s-1) (10-19) (108 yr) (lO^ G) (IQ-i^'s-i) (lO-^^) (10« yr) (10^ G) 



c 


5.757 


-0.91 


< 2.8 


> 3.3 


< 1.3 


1.33 


< 0.24 




> 38 




< 0.4 


D 


5.358 


-0.10 


< 3.0 


> 2.8 


< 1.3 


2.48 


< 1.3 




> 6.7 




< 0.8 


E 


3.536 


2.74 


< 3.0 


> 1.9 


< 1.0 


2.60 


< 1.9 


> 0.05 


> 3.0 


< 110 


< 0.8 


F 


2.624 


2.42 


< 2.1 


> 1.9 


< 0.8 


4.92 


< 1.9 




> 2.2 




< 0.7 


G 


4.040 


-1.09 


< 1.9 


> 3.4 


< 0.9 


4.32 


< 1.3 




> 4.9 




< 0.7 


H 


3.210 


-0.09 


< 1.8 


> 2.8 


< 0.8 


2.21 


< 0.7 




> 7.5 




< 0.5 


I 


3.485 


-1.36 


< 1.5 


> 3.6 


< 0.7 


4.32 


< 1.0 




> 5.3 




< 0.6 


J 


2.101 


-0.51 


< 1.1 


> 3.0 


< 0.5 


1.66 


< 0.24 




> 14 




< 0.2 


L 


4.346 


-2.85 


< 1.3 


> 5.5 


< 0.8 


5.20 


< 1.0 




> 6.8 




< 0.7 


M 


3.677 


-1.08 


< 1.7 


> 3.4 


< 0.8 


1.57 


< 0.18 




> 32 




< 0.3 


N 


3.054 


-0.76 


< 1.5 


> 3.2 


< 0.7 


3.23 


< 0.8 




> 6.4 




< 0.5 


O 


2.643 


1.10 


< 1.8 


> 2.3 


< 0.7 


5.50 


< 1.7 




> 2.4 




< 0.7 


Q 


4.033 


0.80 


< 2.6 


> 2.4 


< 1.0 


1.71 


< 1.0 




> 6.3 




< 0.7 


T 


7.588 


3.84 


< 7.3 


> 1.7 


< 2.4 


4.07 


< 6.0 




> 2.0 




< 2.2 


U 


4.343 


2.15 


< 3.4 


> 2.0 


< 1.2 


1.78 


< 1.7 


> 0.16 


> 4.0 


< 42 


< 0.9 



see Table ^) are all much larger than one would expect for 
similar systems in the Galactic disk, and are therefore, pre- 
sumably, a fossil remnant of gravitational interactions be- 
tween the binary systems and other cluster stars. Rasio fc 



Heggie (1995) calculate the expected value of this eccentric- 



ity for a system with a particular binary period that has 
interacted with other stars in a region of known density for 
a certain length of time. Unfortunately we do not know the 
relevant densities or interaction time scales for the pulsars 
observed in 47 Tuc with any degree of certainty, possibly not 
even to an order of magnitude. Nevertheless, taking plausi- 
ble estimates for density and time scales, we derived some 
"predicted" eccentricities. Those for 47 Tuc E, Q, and T are 
within an order of magnitude of the observed values, which 
is as close to agreement as we can expect to attain given 
the uncertainties in input parameters. The computed eccen- 
tricity for 47 Tuc U is two orders of magnitude below the 
observed value, while that for 47 Tuc H is under-predicted 
by a factor of 1000. 

One possible explanation for the unexpectedly large ec- 
centricity of 47 Tuc U, and possibly that of 47 Tuc H, is that 
they have in the past spent a considerable amount of time 
in a region with a much higher stellar density than they are 
found in at present. For instance, if 47 Tuc L (an isolated 
pulsar) had spent 10^" yr in its present location while part of 
a 2.35-day binary, its computed eccentricity would be 0.15. 
Therefore, some of the pulsars we observe relatively far from 
the centre of the cluster may have non-circular orbits in the 
cluster potential which take them through higher-density 
regions periodically, or may have been ejected from higher- 
density regions to their present locations, leading to some of 
the large eccentricities measured. 

Alternatively, the high eccentricity of the 47 Tuc H bi- 



nary system, by far the largest in the cluster, could indicate 
that it obtained its present companion as an already- formed 
white dwarf through an exchange interaction. However , such 



exchange s tend to produce even higher eccentricities (Phin- 



ney 1992) 



The large eccentricity of 47 Tuc H has permitted a 
measurement of its rate of advance of periastron: u — 
(0.059 ±0.024)° yr"i (Table |). Assuming that this advance 
is entirely due to general relativity, and not to any tidal 
effects (which is hkely, since both stars are presumably de- 
generate and have negligible dimensions compared to the 
orbital separation), the total mass of the binary system is 
1.4lggM0 (all uncertainties here are given at the 2-a level). 
A mass-mass diagram is shown in Fig. |ll|, from which limits 
on the individual masses can be derived: lUp < 2.0 M© and 
rric > 0.1 Mq. It is also improbable that mp < 0.4 M© or 
rric > 0.7 Mq, because this would require cosi > 0.95. 



7 CONCLUSIONS AND PROSPECTS 

In this paper we have presented phase-coherent timing so- 
lutions for 15 millisecond pulsars in 47 Tuc. Our main con- 
clusions can be summarised as follows. 

All pulsars with known timing solutions are located less 
than 1^2 from the cluster centre, and inside this region, their 
spatial density is of the type n{r) cc . This distribution 
is different from that found in M15 {n[r] oc r~'' i, with no 
well-defined outer limit), and possibly similar to the dis- 
tribution of blue stragglers in 47 Tuc with similar masses. 
Their confinement near the centre of the cluster is real, and 
not a selection effect. It is a very important constraint on 
any model of the potential of 47 Tuc. 
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Figure 11. Mass-mass diagram for the 47 Tuc H system. The 
allowed range of total mass (sloping straight line surrounded by 
dashed lines) is derived from the measured rate of advance of pe- 
riastron and its 2-a uncertainty. The area below the thick curving 
line is excluded by the measured mass function and the require- 
ment that cos i > 0. There is a 95% a priori probability that the 
system lies at cosi < 0.95. The vertical lines indicate the average 
measured mass for neutron stars in the radio p ulsar population, 
1.35±O.O4M0 ([Thorsett fc Chakrabarty 1999). 



Two of the pulsars (47 Tuc G and I) have a projected 
separation of only 600 a.u., and similar P/P. Although it 
is unlikely, they may belong to the first system known with 
two detectable pulsars, but further timing measurements are 
required to investigate this. 

We have measured the proper motions of five of the 
pulsars in 47 Tuc. Only the general motion of the cluster 
is detectable, and the value calculated from averaging the 
proper motions of the pulsars is consistent with the proper 
motion of the cluster measured with Hipparcos. 

Lower limits for the surface mass density at several ra- 
dial cuts near the centre of the cluster were derived. We 
have also used the values of (P/P)obs to obtain limits on, or 
estimates of, the accelerations at several projected distances 
from the centre of the cluster. These are consistent with a 
King model of the cluster using accepted values for its dis- 
tance to 47 Tuc, angular core radius and central velocity 
dispersion. 

We inferred upper limits for the magnetic fields of 15 
pulsars: all are below 2.4 x 10^ Gauss. Also, no pulsar has 
a characteristic age smaller than 170 Myr. Additionally we 
determined that the characteristic ages of the pulsars are, 
on average, greater than 1 Gyr. 

The fact that none of the known binary pulsars has an 
orbital period larger than 2.35 days is an important result. 
Sim ulat ions of their formation (Rasio, Pfahl & Rappaport 



200C ) s iggest that we should find binary pulsars with orbital 
periods as large as about one month. However, despite the 



fact that there are no selection effects against the detection 
of such systems, we find none. This suggests that any binary 
pulsars with longer periods were either hardened by interac- 
tions with other stars, disrupted by the high stellar density 
near the centre of the cluster, or were simply outnumbered 
by tight binaries formed at a later stage. The timescale of 
such processes is proportional to the relaxation time, i.e, it is 
proportional to the time between close encounters between 
any two stars. 

Binary formation and hardening is a source of kinetic 
energy for the remaining stars in the cluster and tends t o 
delay (or in some cases reverse) core collapse ( ^pitzer 1987 ); 
the lack of long-period binaries therefore suggests that the 
core of 47 Tuc has had a density similar to or larger than 
the present one for a time that is muc h longer t han the 
relaxation time, i.e., a few billion years (Hut 1992). If the 
cluster had reached the present central density in the last 
few hundred million years, and if it had never experienced 
a high-density phase before, then it should still contain bi- 
nary pulsars with long orbital periods. As we have also seen, 
the ages of the pulsars are of the order of a few Gyr, and 
their large numbers in this cluster suggest that vigorous mil- 
lisecond pulsar formation occurred a long time ago, which 
is consistent with the idea that the core of 47 Tuc was once, 
or has been for a long time, in a high-density state. 

Another piece of evidence of possible relevance to this 
question comes from the higher than expected eccentricities 
of binaries like 47 Tuc U and H, the latter of which enabled 
us to measure its rate of advance of periastron, yielding a 
total mass rup + rric — 1.4lgg Mq. These eccentricities sug- 
gest that the product of stellar density and velocity disper- 
sion near these binaries was higher than is observed today 
for a long time, possibly several Gyr (§ Whether this was 
so because the pulsars spent a considerable amount of time 
nearer the present-density core, or because the core was once 
denser, it would imply a dense core for a very long period 
of time. 

If the core of 47 Tuc was in fact relatively dense for a 
few Gyr, it is likely that this contributed to the large num- 
bers of millisecond pulsars in the cluster (via the mecha- 
nisms of binary formation and hardening in dense cores) . In 
that case the number of millisecond pulsars in a cluster de- 
pends heavily on its previous dynamical history. This could 
explain why clusters with similar properties (core density, 
total mass, stellar mass distribution, metallicity) have large 
differ ences in their pulsar populat ions, in both number and 
kind ( Kulkarni fc Anderson 1996 ). More complete searches 
for pulsars in globular clusters, a better understanding of 
selection effects and a good determination of the evolution- 
ary states of globular clusters will be needed to address this 
question. 

A continuing timing programme will lead to additional 
and improved measurements of the proper motions of the 
pulsars. A measurement of the proper motions of the pulsars 
relative to the cluster would add constraints to the cluster 
mass model. Continued timing will also allow a better de- 
termination of the orbital parameters of 47 Tuc H, and will 
test whether 47 Tuc G and I are physically associated. Addi- 
tional solutions for known pulsars may also be determined, 
and new pulsars will almost certainly be discovered. Contin- 
ued monitoring of this cluster is therefore likely to remain a 
worthy enterprise. 
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